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Abstract 

We consider covariant metric theories of coupled gravity-matter systems satisfying the fol- 
lowing two conditions: First, it is assumed that, by a hyperbolic reduction process, a system of 
first order symmetric hyperbolic partial differential equations can be deduced from the matter 
field equations. Second, gravity is supposed to be coupled to the matter fields by requiring 
that the Ricci tensor is a smooth function of the basic matter field variables and the metric. It 
is shown then that the "time" evolution of these type of gravity-matter systems preserves the 
symmetries of initial data specifications. 

1 Introduction 

In a recent paper ||l9| the necessary and sufficient conditions ensuring the existence of Kilhng vector 
fields in covariant metric theories of coupled gravity matter systems were investigated. In that work 
matter fields were assumed to satisfy quasi-linear wave equations and to possess minimal coupling 
to gravity. In addition, gravity was supposed to be coupled to the matter fields by requiring that the 
Ricci tensor is a function of the basic matter field variables, their first covariant derivatives and the 
metric. Within this framework it was demonstrated that symmetries of initial data specifications 
are preserved by "time" evolution. 

The pragmatic aim in jl^ was to justify some of the claims of |^, ^ concerning the existence of 
Killing vector fields in the characteristic initial problem in case of Einstein-Klein-Gordon, Einstein- 
[non-Abelian] Higgs or Einstein- [Maxwell]- Yang-Mills-dilaton systems. The matter field variables^] 
of these particular Einstein-matter systems all satisfy quasi-linear wave equations. However, the evo- 
lution of various matter fields, such as perfect- or dissipative fiuids and spinor fields, is not governed 
by quasi-Hnear wave equations; rather they satisfy first order symmetric hyperboHc equations. 

Clearly, it is of obvious physical interest to enlarge the framework of to include matter fields 
satisfying equations of this more general type and to determine those conditions which guarantee 
the existence of a Killing vector field for the associated gravity-matter systems. This is, in fact, 
our principal aim in this paper which is organized as follows. The next section is to specify the 

*This research was supported in part by OTKA grant T034337. 

^In case of gauge fields they are yielded by suitable gauge fixing (see e.g. [hs, hsj for further details). 
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class of gravity-matter systems to which our main result applies. In particular, a representation of 
the physical fields and the relevant basic field equations are chosen so that well-posed initial value 
problem can be associated with the selected gravity-matter systems. Section || is to study all the 
equations relevant for the 'evolution' of the Lie derivatives, with respect to a 'candidate KiUing 
vector field', of our basic variables. This section includes our main result justifying the claim that 
the "time" evolution of the selected gravity-matter systems preserves the symmetries of initial data 
specifications. Finally, section ^ contains our concluding remarks. In particular, issues associated 
with the evolution of matter fields on fixed geometrical backgrounds and the possible affection of a 
non-dynamical field on the evolution of symmetries are considered. 

2 The gravity-matter systems 

In order to have a suitable framework for the study of spacetime symmetries, first of all we give 
a mathematically precise specification of the gravity-matter theories to which our main result will 
apply. Obviously, it is impossible to recall details of all the particular systems which have ever been 
studied. Instead, a general, and thereby fiexible, enough framework is chosen within which common 
features of a reasonably large set of particular theories can be investigated. The model specified 
below is highly infiuenced by the fundamental study of Geroch [0, where a universal treatment of 
physical fields and equations is presented. 

2.1 The choice of basic variables and field equations 

The physical fields are supposed to be represented by smooth tensor fields on a smooth four- 
dimensional paracompact orientable manifold M. In particular, the "gravitational field" is assumed 
to be represented by a smooth Lorentzian metric field gab on M and all the gravitational effects 
are expected to be described in terms of this metric and covariant quantities derivable from it. It 
is known that other basic variables such as tetrads or spin frames could also be used to represent 
gravity. However, since our eventual aim in this work is to study spacetime symmetries the use of 
gab suits the most.0 

In general, tensor fields representing the matter content of a theory will be denoted by T^^, 
where i in round brackets is a name index taking values from the set {1,2, ...,7V}, and the capital 
script letter A stands for a composite (abstract) index. The constituents of A, i.e. the exact form 
T,.)' " '''j^^ of the {ki,li) type tensor field T^-^ will be spelled out explicitly in all the non-self- 
explaining situations. The matter fields might also have gauge dependence but, even if they have, it 
will be assumed that suitable gauge choices have been made so that they can be considered as fields 
on the spacetime manifold. Thereby the relevant gauge or internal space indices will be suppressed. 

The evolution of the fields T-^ is supposed to be given in terms of partial differential equations. In 
particular, T"^ are assumed to solve systems of first order quasi-linear partial differential equations 
of the form 

E<;,«v.r^'+^;!,=o, (2.1) 

were Va denotes the unique (torsion free) covariant derivati ve o perator associated with the spacetime 
metric gab- Moreover, the coefficients K'^.^^,^^ and in ( |2.l| ) are expected to be smooth fields on 
a tensor bundle, T, built up from the fields T-^ and gab- 

^Unless otherwise stated we shall use the notation and conventions of [El|. 
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It is worth emphasizing that this framework covers matter systems Uke Klein-Gordon fields, 
Maxwell fields, perfect fiuids and various other type of matter fields, along with their conventional 
couplings (see for more details Appendix A of ||ll|). In fact, it seems to cover virtually every partial 
differential equation in physics.^ In particular, quasi-linear wave equations of the form 



V/T« + (g,,, T^^ , V.T^^ ) = 0, (2.2) 



can immediately be put into the form of (|2.l|), by following a standard procedure (see e.g. pages 
76-77 of iQ) which starts, by adding the first derivatives T-^^.^e — VgT"^ to the basic variables. 

Note also that, as indicated above, in case of matter fields, like Yang-Mills or Higgs fields with 
a gauge freedom, to get the desired framework a gauge fixing is required to be made. For instance, 
in case of a Yang-Mills field the use of the generalized Lorentz gauge condition is needed to put 
the Yang-Mills equations into the form of a quasi-linear wave equation (see for details e.g. p^), 
which can then be recast into the form of ( |2.l[ ) by the above mentioned procedure. 

According to the above assumptions matter fields are assumed to be coupled to gravity through 
the dependence of K-^^^ _^ b and on the metric, as well as, via the presence of the metric compatible 



covariant derivative operator Vg in (p.!]). Concerning the coupling of the geometry to matter fields, 
we shall require that the Ricci tensor Rab be a smooth function of the matter field variables and the 
metric, 

Rab = Rab{Tf', gmn) ■ (2.3) 



Note that this assumption immediately imposes a restriction, in addition to (2J), on the fields 
via the twice contracted Bianchi identity. 

It is straightforward to see that the conditions so far have been made are satisfied by all the 
'customary' Einstein-matter systems. In particular, the present framework immediately covers the 
setting applied in |[l9j since the matter field variables used here includes the first covariant derivatives 
of the basic matter field variables applied there. Moreover, the selected model is general enough to 
host e.g. the 'conformally equivalent representation' of higher curvature theories associated with a 
gravitational Lagrangian that is a polynomial of the Ricci scalar (for further details see e.g. [p^). 

In addition to the above assumptions we shall impose the following requirement which plays a 
significant role in the derivation of our main result and which is usually made implicitly in most of 
the related discussions. 

Condition 2.1 The fields K-^'^ n and L-^ , alonq with the Ricci tensor when it is viewed as a 



functional of the basic variables as the expression appearing on the r.h.s. of p. j| A are assumed to 
be smooth functionals on the tensor bundle T depending exclusively on T-^^ and gab, 



Ae Ti'Ae 



B{Tf^ygmn) and Lf^^=L^^^{Tf^^,gran), (2.4) 



without having any explicit dependence on the points of Ad . 
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In principal, spinor fields could also be present amonty the allowed matter field variables. This is particularly 



because the Dirac equations possess the form of (2.1). Moreover, it was demonstrated by Priedrich and Rendall 
that symmetric hyperbolic equations can be built up from the Einstein-Dirac system. The main purpose of this work, 
however, is to study the existence of symmetries which requires the use of Lie derivatives of all the fields with respect 
to essentially arbitrary vector fields on M . On the other hand, it is known that the Lie derivative of spinor fields 
cannot be defined (at least not in the conventional way) with respect to arbitrary vector fields but only to conformal 
Killing vector fields. Before the procedure described below can be extended to Einstein-Dirac systems as well one 
should clear up e.g. how to make a sensible variation of a composite metric-spinor system. Obviously, there is no 
problem of varying the spinor field while the metric is kept fixed. However, a procedure of complete generality, such 
as 'changing the metric while keeping the spinor field fixed' does not exist. 
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This assumption expresses the idea that the distinction between two points of a spacetime man- 
ifold M is rooted in the difference between the values of some physical field at that points. In other 
words, the points of the spacetime {M,gab) are assumed to be identified by observing some physical 
fields not by somehow "perceiving the points of M themselves directly" ||ll|. Notice, however, that 
whenever the dynamics of certain fields T^^ on a fixed background spacetime is considered the co- 
efficients K-^^^_^is and will immediately pick up an 'explicit' spacetime dependence through that 
of gab- 



Note finally that the relations (2J.) and (p.3|), along with condition 2A, manifest an important 
property of the composite gravity-matter system. Namely, they guarantee that the coupHng of all 
the subsystems rests upon only zero order terms of the relevant basic variables. 



2.2 Hyperbolic reduction of the field equations 

A key feature of the field equations of physics is that a system of partial differential equations can 
be deduced from them so that, for suitable initial data specifications, this equations yield the 'time' 
evolution of the relevant system. This subsection is to impose further conditions which guarantee 
the well-posed initial value problems can be associated with the selected gravity-matter theories. 



2.2.1 Symmetric hyperbolic systems for the matter fields 



To attain 'existence and uniqueness results' relevant for an initial value problem associated with the 
gravity-matter systems of the above type first the basic field equations have to be put into a specific 
form of PDEs for which existence and uniqueness of solutions is guaranteed. In this work first order 
symmetric hyperbolic PDEs are used for this purpose which are specified as follows |p^: 

Consider, as our unknown, a function $ defined on an open subset of R** and assume that it takes 
values in a finite dimensional vector space V . The hyperbolic systems we are interested in now are 
of the form 

A^d-y<^ + B = 0, (2.5) 

where A"! — $) are supposed to be smooth functions defined on an open subset of x F 

taking values in the vector space L{V) of Hnear maps of V to itself, while B — B{x, $) is expected 
to be a smooth function on an open subset of R'* x F with values in V. The basic requirement on 
the system to be called symmetric hyperbolic is that the A'^ are symmetric, with respect to some 
inner product on V, and there exists a covector field ^-y such that A^^t is positive definite for all 
admissible values of <&. 

By imposing the following condition we shall restrict our considerations to matter fields for which 
a symmetric hyperbolic system can be deduced from 



Condition 2.2 Equation (2.1) can be split into a symmetric hyperbolic system of the form (2.5), 
in terms of a suitable variable $ labeling cross-sections of the tensor bundle built up from the fields 
T^^, and constraint equations which, along with the possible gauge conditions, propagate under the 
evolution of the associated hyperbolic system. 



Note that there is no known a general practical procedure which would immediately produce a 
symmetric hyperboHc system of the form (2.5) from (2.1). In fact, the actual way of deducing a 



symmetric hyperbolic system (2.5) from ( |2.1j ) and the justification that the constraints also propagate 
for each particular case might be very complicated and most of the time requires very special 
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considerations. However, the actual way of getting the required hyperbolic reduction of the selected 
matter system will not play any role in the basic argument concerning the evolution of symmetries 
of initial data specifications. 



Notice also that even after imposing condition 2.2 the framework we have still covers matter 
systems like Klein-Gordon, dilaton, Higgs, Maxwell- Yang-Mills fields fields and various fiuids along 
with the conventional couplings of these systems. In fact, this setting is suitable to host virtually 
all the physically relevant matter systems for which a well-posed initial value problem can ever be 
formulated. 



2.2.2 Hyperbolization of the gravitational part 



To get symmetric hyperbolic equations relevant for the selected gravity-matter system we need to 
take care beside the hyperbolic reduction of (^) that of (2.3), as well. It is known |11, 1^ that, due 
to the general covariance of the considered metric theories of gravity, there is no way to get hyperbolic 



evolution from the coupled field equations (2.1) and (p7S) without fixing the gauge associated with 



the diffeomorphism invariance. One of the possibilitiesr 
form of the Ricci tensor Rab given in local coordinates x 

1 



to make such a gauge fixing rests on the 

' as 



Rati = "l^g^^'dpO^gaP + gS(adp)T + Ha[j [gep, d^Qep) , 



(2.6) 



and it is done by simply replacing T° = g^fV^p in ( |2.6D with arbitrary 'gauge source functions' / 
(chosen so that = is satisfied on initial hypersurfaces) . The equations yielded by substituting, 
in local coordinates, the r.h.s. of (2.6) for the l.h.s of (|2.3|), with replaced by f^, are refereed to 
as reduced gravity equations. These equations can be put into the form of a first order symmetric 
hyperbolic equation by introducing all first order derivatives of ga/3 as new variables and let $ consists 
of 5a/3 together with these derivatives. Then by making use of the reduced gravity equations a system 
of symmetric hyperbolic equations for <i> can be built up by a standard procedure (see for more details 
e.g. illl). 



2.2.3 Hyperbolic reduction of the coupled system 

Upon having a hyperbolic reduction of the matter and gravity equations separately, it is important 
to know whether the coupled equations also form a symmetric hyperbolic system. The answer is 
in the affirmative. This can be justified by recalling that the coupling of the matter field equations 



(2.1) - which by condition 2.2 can be cast into a symmetric hyperbolic system of the form (2.5) 



and that of a fixed hyperbolic reduction of (|2.3D - written in first order symmetric hyperbolic form - 
is only by terms of order zero, hence, the combined system is also symmetric hyperbolic. Therefore, 
local existence and uniqueness results immediately apply to the relevant reduced gravity-matter 
equations in case of regular initial value problems. 

There is only one significant requirement the applied initial value problem has to satisfy: The 
existence and uniqueness of solutions to symmetric hyperbolic equations is expected to be guaranteed 
within its framework in the smooth setting. The initial value problems which satisfy this condition 
are referred to as regular initial value problems. Immediate examples for appropriate initial value 
problems are the standard Cauchy problem (see e.g. Ref. |Q|) and also the characteristic initial 
value problem associated with an initial hypersurface represented by either the union of two smooth 
null hypersurfaces intersecting on a 2-dimensional spacelike surface 0, Q or a characteristic cone 

''For an excellent review of all the related issues see jsj]. 
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1^, |l|, ^ . Since no further requirement on the initial value problem is used anywhere in the derivation 
of our results we shall not make a definite choice among these regular initial value problems. Thereby 
the following notation will be applied on equal footing to any of these initial value problems: The 
initial hypersurface will be denoted by S, while the initial data on E will be represented by the basic 
varia bles in square brackets. In particular, an initial data specification for an equation of the form 
(2^) is represented by the pair (S, [$]). 

As mentioned above the hyperbolic reduction is done by splitting the basic equations to pure 
evolution equations and constraint ones. In fact, the constraint equations, along with the gauge 
conditions, manifest certain (differential) identities that have to be satisfied by T^^ and gab- In 
particular, they impose restrictions on the possible initial data sets. Thereby, having a hyperbolic 
reduction of the coupled system suitable initial data consistent with the relevant gauge conditions 
and satisfying all the constraints has to be arranged. Finally, it has also to be demonstrated that 
the gauge conditions and the constraints propagate with the 'time' evolution of the relevant coupled 
reduced system. In the present case, by making use of a suitable combination of condition 2^ and 
the twice contracted Bianchi identity (see also page 20 and pages 76-77 of it can be justified 

that the gauge conditions and the constraints will be satisfied throughout the domain of dependence 
associated with a solution to the coupled evolution equations if they hold on the initial data surface. 



3 The evolution of symmetric initial data specifications 

Consider a gravity-matter system of the type specified in the previous section satisfying conditions 



2.1 and 2.2. Assume that a regular initial value problem has been chosen, moreover, a suitable 
initial data set relevant for the coupled symmetric hyperbolic evolution equations also satisfying the 
pertinent gauge conditions and constraint equations has also been selected. Recall that then we 
immediately have a detailed knowledge of the metric and matter fields on the initial data surface 
S since in any of the regular initial value problems T^-^ and gab, along with their derivatives up 
to arbitrary order, can always be determined on E by making use of an initial data specification, 
(S, [T^^], [gab]), and the field equations. Clearly, it is of obvious interest to identify those initial data 
sets that are guaranteed to evolve to configurations possessing certain type of symmetries. The rest 
of this paper is to make it clear under what conditions the symmetries, and more importantly what 
sort of symmetries, of initial data sets will be inherited by the corresponding solutions. 

There are various possibilities which make it sensible to consider an initial data specification to 
be symmetric. The most straightforward one is the following: An initial data set (E, [T^^], [gat]) is 
considered to be symmetric if there exists a diffeomorphism i/i : S ^ E of the initial hypersurface 
S onto itself so that it leaves the initial data invariant. These type of symmetries of initial data 
specifications were considered in jl^, for the case of Einstein-vacuum systems. Based essentially on 
the uniqueness of maximal Cauchy developments, it was demonstrated there that such a symmetry 
of a vacuum initial data set will be inherited by the corresponding solution. By making use of a 
straightforward adaptation of this argument (see, for more details, pages 90-91 of |Q) to the case of 
gravity-matter systems investigated here, it can also be justified that a symmetry, of the considered 
type, of an initial data set gives also rise to a symmetry of the associated Cauchy development, i.e. 
there will exist a diffeomorphism : M M with = and ^'*5ah = gab with ^fjs = ip- 

Since there are symmetries of Cauchy developments which does not necessarily map any initial 
data surfaces onto themselves it is of obvious interest to know whether the traces of these more 
general type of symmetries can also be recognized by making use of initial data specifications. It 
seems to be hopeless to find an evidence for the existence of these more general type of symmetries in 
the case of discrete symmetry transformations. However, whenever a gravity-matter system admits 
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a Killing vector field, i.e. there exists, at least locally, a one-parameter family of isometry actions in 
a neighbourhood of S, the procedure described in detail in the following section can be applied. 

Before approaching to this point there are two comments in order. Firstly, I would like to 
emphasize the complementary character of the procedures associated with the above mentioned two 
types of symmetries. The results relevant for the more general type of symmetries obviously cover the 
special case whenever an initial data specification admits a one-parameter family of diffeomorphisms 
-i/'r : S ^ S leaving the initial data invariant. However, it has its drawback not being able to handle 
the case of discrete isometry actions which, on the other hand, can be covered by the indicated 
generalization of the method of Priedrich and Kendall (l^ in the particular case when the initial 
data surface is mapped onto itself by the associated diffeomorphism. 

Secondly, the results presented below are due to straightforward adaptations of earlier methods 
applied in the Einstein- vacuum case to study the linearization stability problem by Moncrief [l^ 
(see also [|[ §])■ Note that these results provide an immediate generalization of the main assertion of 
p9| , and motivates a possible further strengthening of the conclusion of jTs), since the results covered 
by ||l9t [T^ are relevant merely for gravity-matter systems with matter fields satisfying quasi-linear 
wave equations. 



3.1 The construction of a candidate Killing vector field 

It is known that if we had a Killing vector field if", satisfying the Killing equation, 

C-KQab = VaiiTb + VbKa = (3.1) 

and the integrability condition, 

VaVfeX, + RbcjKd - 0, (3.2) 

of the 2-form field Eab — ^ aKb this Killing vector field would be completely determined by the 
values of Ka and VaKb at any point of M since the above two equations imply a system of ODE's 
for the components of Ka and VaKb along any curve. However, the existence of a Killing vector 
field cannot be proven in this way. Nevertheless, the contraction, 

V^VeiC" + R'^dK'^ = 0, (3.3) 



of (3.2), which is a linear homogeneous wave equation for if, provides the means to construct a 
'candidate' Killing vector field. Obviously, any Killing vector field satisfies (3.3) but not all of its 
solutions will give rise to a Killing vector field. 



3.2 Evolution of the Lie derivatives of the basic variables 

The conditions ensuring the existence of initial data [K°-] on an initial hypersurface S so that the 



unique solution of (3^) will be a Killing vector field are necessarily given in terms of restrictions 
on the values of the fields T^^ and gab on E. The relevant requirements can be read off the evolution 
equations for the Lie derivatives of our basic variables. The remaining part of this section is to 
identify these conditions and thereby we shall prove the following: 

Theorem: Let {M, gab) be a spacetime associated with a gravity-matter system as it was specified 
in section^. Denote by D[S\ the domain of dependence of an initial hypersurface E associated with 
a regular initial value problem. Then there exists a non-trivial Killing vector field K^ on -D[E], so 
that the matter fields are also invariant, i.e. CkT^^ = 0, if there exists a non-trivial initial data set 

for ^ so that jOkT^^, ^Kgab and VciCKgab) all vanish identically on E. 
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Before presenting the proof of this theorem let us emphasize certain points relevant for the 
applied setting. First of all, the practical problem of finding initial data [K""] for ( ^7^ ) - which in 
each particular case is the hard part of the work - is not considered here. It is, however, important 
to note that, upon having any sort of initial data [K'^] for (^), the validity of our conditions, i.e. 
the vanishing of CkT-^ CKQab and VciCxgab), can be immediately justified. This is a consequence 



of the fact that knowledge of [K°-], along with the equation ( |3.3D , suffices to determine derivatives 
of K'^ up to any order on E. As already mentioned above, the fields T"^ and gab, along with their 
derivatives up to arbitrary order, can also be determined on the initial hypersurface S in any regular 
initial value problem. Thereby the existence of a local one-parameter group of diffeomorphisms 
acting on the Cauchy development which leaves the matter fields and the geometry invariant can be 
read off the initial data for T^^ and gab in advance of having the relevant solution explicitly. 

Secondly, in case of gauge fields it might seem to be too strong to require the vanishing of the Lie 
derivative of the field variable. It is known, for instance, that in the particular case of a Yang-Mills 
field the invariance of a vector potential Aa (taking values in a Lie algebra q) under the action of a 
one-parameter group of diffeomorphisms associated with a vector field if implies only that CxAa 
is equal to the gauge covariant derivative Da of a g-valued function W, i.e. CxAa = DaW. This 
general invariance property plays an important role whenever there are several symmetries acting 
simultaneously on a given spacetime. Note, however, th at in the present situation - when we look 
for a single candidate Killing vector field determined by (3^) and suitable initial data specifications 
- the function W can be gauged away (see, e.g., page 19 of 1 1|) so that the invariance of an adapted 
vector potential Aa can, in fact, without loss of generality, be expressed as CxAa = 0. 



Proof: To start off consider a vector field K'^ satisfying ( |3.3| ) but which is kept otherwise to be 
arbitrary. Then, by taking the covariant derivative of ( |3.3|) , commuting derivatives and applying the 
contracted Bianchi identity, it can be shown that Cxgab satisfies the equation 

2R\b^iCKgef) 



V^Ve {Cxgab) — —^CxRab 

Moreover, by taking the Lie derivative of (2^) we get^ 



2R\^{CKgb)e). 



(3.4) 



^KRab 




dRat 
dg. 



i^Kg. 



Then, in virtue of ( p. 4 ) and ( |3.5|) , Cxgab satisfies an equation of the form 

V^Ve {CKgab) = Pab{CKgcd) + Q i., abiCxT;^^) 

where Pab and Q,.,afc are linear and homogeneous functions of their indicated arguments. 



(3.5) 



(3.6) 



The next step is to show that CkT^^ satisfy an equation analogous to (2.1) so that the coefficient 
of the principal part of the relevant equation is equal to that of (2J_). To achieve this take the Lie 
derivative of (2J_) with respect to the vector field if". This yields the relation 




dg, 



cd 



^Kg 



cd 



> CKgcd = 0. 



(3.7) 



'whenever T-^ and S'^ are tensor fields of type (k,l) and (m, n), respectively, (dT^ /dS^) is considered to be a 
tensor field of type {k + n,l + m). Accordingly, the contraction {dT-^ /dS'^)CK S'^ is again a tensor field of type {k,l). 
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To put ( p?^ ) into the form of (2.1) the commutation relation of the operators Ck and Va is needed 
to be used which reads as 



Ck (V,r^) = Ve [CkT^) + £ (rt!) U [V^i^ff] :\ - E (r-:;) [-^CKg] C (3-8) 

where (r-^^) fa,] and (t^) i^*' stand for r,,)"'"2"°""f,i...6,^ and T(„ I. .6,^ , respectively; 
moreover, the notation 



(3.9) 



has been applied. Then, in virtue of ( p.Tj ) and ( p. 8]) , the desired equation reads as 

E<:.-^e(/:Kr«)+L;^,=o, (3.10) 



0) 



where the coefficients of the terms Ve(£i<-T^^j) are exactly that of VgT^j in (2J.), and, the stand 
for a sum of terms which are linear and homogeneous expressions of either of the variables LkT'^^ , 
CKQab or VciC-Kgab), respectively. 

To complete our argument notice first that, in virtue of the coincidence of the principal parts of 
(3.10) and ( |2.1[ ), along with condition 2.2, the same hyperbolization procedure that was supposed to 
exist in case of (EJ) can also be used to deduce from (3.10) a symmetric hyperbolic system for the 



variables £i<-T-^ rlSimilarly, the linear wave equation (3.6) can be put into the form of a symmetric 



hyperbolic system for the variables Cxgab and VcCxgab- Since the coupHng of these two systems 
happens only through terms of zero order the coupled system immediately gives rise to a symmetric 
hyperbolic system for the fields CkT^^, CRgat and '^dCKgat)- In addition the 'source term' of 
the resulting system is a sum of expressions which are linear and homogeneous in either of these 
variables. This property guarantees that the combined equation necessarily have the identi cally zero 
s oluti on for vanishing initial data. Since the symmetric first order system derived from (3^) and 
( 3.10 ) is Hnear and homogeneous in the basic variables the constraints, that there might be yielded 
by the relevant hyperboHc reduction of these equations, are immediately satisfied by the identically 
zero solutions of the coupled evolution equations. Hence whenever initial data [if] can be chosen 
on an initial hypersurface S so that CkT^^, C-Kgab and '^dC-Kgab) are zero on S then each of the 
fields £kT^^ and Lxgab will identically vanish throughout the domain where the associated unique 
solution of (3^) does exist. To see, finally, that this domain has to coincide with D[E] note that, 
since ( |3.3| ) is a linear homogeneous wave equation any solution of it can be shown - by making use 
of the 'patching together local solutions' techniques described e.g. on page 266 of [|6) - to extend 
over the entire domain of dependence Z?[S] of S. q 



4 Concluding remarks 



In virtue of the above result one may conclude that for any gravity-matter systems satisfying condi- 
tions 2A and 2^ Killing symmetries of initial data specifications (S, [T^^], [gab]) are preserved by the 



^Recall that a hyperbolic reduction of an equation of the type (2.1), whenever it can be made, is achieved by 
a regular algebraic procedure (see e.g. page 9 of |llj|) which depends on the principal part of the equation exclu- 
sively. Thereby, however different the source terms and might be, their difference has no influence on the 



hyperbolization of equation (3.10) 
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"time" evolution. This section is to consider situations where either of the conditions of the above 
theorem fails to be satisfied. 



Let us start by considering the evolution of certain fields on a fixed non-dynamical spacetime 



{M,gab)- Suppose that if is a vector field on M satisfying (3.3) but the Lie derivative of the metric 
with respect to if should not necessarily vanish. It is straightforward to see that whenever CxQab 
does not vanish the relevant form of (B^) gives rise to a constraint equation which has to be satisfied 



by the Lie derivatives of fields T-^^. Furthermore, (3.10) determining the evolution of £kT^^ will 



have non-trivial source terms in consequence of the presence of the non-vanishing fields Cxgcd and 
VeC-Kgcd- This impHes then that despite of the fact that the fields T^^ might be arranged to possess 
certain type of symmetries on the initial hypersurface the evolution of the system will suppress this 
property unless if is a Killing vector field of the fixed background spacetime. 

The situation is very similar, although the equations simplify considerably, whenever if is chosen 
to be a Killing vector field for the fixed background geometry but the matter fields are kept to be 
general. Then the non-invariance of the matter fields - represented by the Lie derivatives ^kT^^ - 
'evolve' according to a linear homogeneous equation of the form 

E^^fx^^V. +^i.- , (r.T^,) ^ 0, (4.1) 

(J) (fc) 

where L-f^^ts are some smooth fields, moreover, the Lie derivatives CkT;^ have also to satisfy the 

{k) (ij 

additional constraint 



'dRab\ 



CkT^^ = 0. (4.2) 



An analogous argument applies whenever the geometry is chosen to be dynamical and initial 
data is arranged for (3^) so that CxQcd and VeC-KQcd vanish identically on E, however, either of 
the Lie derivatives £kT^^ is not zero there. Then the unique solution K'^ of (3^) will not give rise 
to a Killing vector field of the coupled gravity-matter system and, in general, neither of the fields 
CkT-^ , LxQcd and VeCxgcd will vanish on succeeding Cauchy surfaces. 



Consider now a situati on whenever condition 2.1 is not satisfied, i.e. suppose that either of the 
fields K-^^B, L-^ used in (2J) or the expression appearing on the r.h.s. of (2^) depends expHcitly 
on the points of the spacetime. This can occur, for instance, whenever there is an additional field 
present not taking part in the evolution. It is straightforward to see that then the Lie derivative of 
either of the fields K-^^^, L-^ or Rab will not be a homogeneous expression of the Lie derivatives of 
the dynamical fields exclusively. Instead, terms containing the Lie derivative of the non-dynamical 
field appear in equations (3.6) and ( B.1C| ). Hence, unless this additional non-dynamical field is itself 
invariant under the action of the one-parameter famihes of diffeomorphisms associated with our to 
be KiUing field if, these terms behave as source terms yielding the loss of any possible KiUing 
symmetry of an initial data specification. 

It is important to emphasize that the presented results also provide means to justify or dispel 
certain assumptions (usually made implicitly) which are applied e.g. in searching for exact or 
numerical solutions of the type of gravity-matter systems investigated here. For instance, if a 
system does not satisfy condition |2.l| it straightforwardly follows from the above considerations that 
the application of a dimensional reduction, based on the assumption that the spacetime admits a 
global symmetry, might not be completely consistent with the system under investigation. Clearly, it 
may happen that a gravity-matter system of this type has no solution which would possess a (local) 
one-parameter group of isometry actions which also would leave the matter fields to be invariant. 
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Let us finally mention that on a similar base the results presented in this paper might be useful in 
testing the trustworthiness of a code in numerical simulation of gravity-matter systems. In particular, 



in case of theories satisfying conditions 2.1 and |2.2| one can be sure that something has to be wrong 



with a code if it does not preserve (at least approximately) symmetries of initial data specifications. 
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